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Abstract 

_ , In Montgomery (2009a), we show that accretion discs in binary systems could retrogradely precess by 

' the same physics that causes the Earth to retrogradely precess [i.e., tidal torques by the Moon and the 

Sun (or secondary star) on a tilted, spinning, non-spherical Earth (or e.g., primary star surrounded by an 
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accretion disc)]. In addition, we show that the state of matter and the geometrical shape of the celestial 
object could significantly affect the precessional value. For example, a Cataclysmic Variable (CV) Dwarf 
Novae (DN) non-magnetic system that shows negative superhumps in its light curve can be described 
lO \ by a retrogradely precessing, difi'erentially rotating, tilted disc. Because the disc is a fluid and because 

the gas stream overflows the tilted disc and particles can migrate into inner disc annuli, coupled to the 
disc could be a retrogradely precessing inner ring that is located near the innermost annuli of the disc. 
' However, numerical simulations by Bisikalo et al. (2003, 2004) and this work show that an inner spiral 

\ density wave can be generated instead of an inner ring. Therefore, we show that retrograde precession in 

non-magnetic, spinning, tilted CV DN systems can equally be described by a retrogradely precessing and 
differentially rotating disc with an attached retrogradely precessing inner spiral density wave so long as 
the wave appears at the same radius as the ring and within the plane of the tilted disc. We find that the 
theoretical results generated in this work agree well with the theoretical results presented in Montgomery 
(2009a) and thus with the numerical simulations and select CV DN systems in Montgomery (2009b) that 
^ \ may have a main sequence secondary. Therefore, pressure effects do need to be considered in CV DN 

systems that exhibit negative superhumps if the accretion discs are tilted and have an inner spiral density 
wave that is in the plane of the disc. 

\^ Keywords: accretion, accretion discs - binaries: general - binaries: close - stars: dwarf novae - novae, 
: cataclysmic variables - hydrodynamics 
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1 Introduction 



Variability in non-magnetic. Cataclysmic Variable (CV) close binaries has been largely attributed to in- 
stabilities. The outbursts in CV recurrent novae (RN) and CV classical novae (CN) can be caused by 
thermonuclear runaway on the surface of a white dwarf or by a thermal and viscous instability in the disc. 
^ ' The latter type of outburst involves cycling from low-to-high thermal and viscous states (see e.g., Meyer & 
. Meyer-Hofmeister 1981 and references within), and thus cycling from low-to- high mass transfer rates (see 
e.g., Lasota 2001), that results in a release of accretion energy as mass is drained onto the white dwarf. 
Transient X-ray binaries with neutron stars or black hole primaries and CV Dwarf Novae (DN) are also 
subject to thermal and viscous instabilities. However, the thermal and viscous instabilities in CV RN are 
larger in scale than those in CV DN. In addition to thermal and viscous instabilities, CV DN can experience 
a tidal instability. Results are CV DN SU UMa and Nova-Like (NL) outbursts. The larger amplitude and 
longer duration superoutbursts involve both a thermal-viscous instability and a tidal instability. 

Both CV DN SU UMa and NL outbursts and superoutbursts show characteristic hump-shaped modula- 
tions, known as positive superhumps, in their light curves. These positive superhump modulations have a 
period that is a few percent longer than the orbital period. They could be caused by dissipation from spiral 
shocks that formed from converging gas flows in the middle and outer annuli of the disc (see e.g., Smith et 
al. 2007). The source to this spiral shock generation is the j : j — 1 = 3 : 2 resonance (Lubow 1991a, b), using 
the terminology established by Frank, King, & Raine (1992). A competing theory has also been suggested 
by Bisikalo et al. (2004) - the superoutburst is associated with the generation of an spiral density wave in 
inner annuli of the disc, a zone that is denser and gas-dynamically un-perturbed. This inner spiral density 
wave contributes to a substantial increase in mass transfer rate in that region. 



In addition to positive super humps, some non- magnetic CV DN systems show negative super humps 
in their hght curves. These negative superhump modulations have a period that is a few percent shorter 
than the orbital period. Some systems' light curves only show negative superhumps (e.g., CN Ori), or 
only show positive superhumps (e.g., SU UMa), or show both modulations simultaneously [e.g., SDSS 
J210014. 12+00446.0 (Olech et al. 2009)], or show one type of modulation only later to be replaced by 
the other [e.g., TT Ari (Kim et al. 2009)]. Negative superhumps are suggested to be a consequence of a 
partially tilted disc (Patterson et al. 1993) or arc related to a warped accretion disc (Petterson 1977; Murray 
& Armitage 1998; Terquem k Papaloizou 2000; Murray ct al. 2002; Foulkes, Haswell, & Murray 2006). Our 
previous numerical simulation work (Wood, Montgomery, &; Simpson 2000; Montgomery 2004, Montgomery 
2009b) shows that negative superhumps can be produced by a fully tilted disc, and these modulations have 
a shape and a period that are similar to those in observational light curves. 

Accretion discs can warp or tilt via a potpourri of suggested sources. For example, disc tilt in X-ray 
binaries can be from gas streaming at an upward angle from the inner Lagrange point for one half of the 
orbit and at a downward angle for the second half of the orbit (Boynton et al. 1980). In some CVs, a 
disc tilt can be held constant by gas stream being fed by the magnetic field of the secondary (Barrett, 
O'Donoghue, & Warner, 1988). Also for CV systems, a disc tilt instability can result from a coupling of 
an eccentric instability to Lindblad resonances (Lubow 1992). A vertical resonant oscillation of the disc 
mid-plane can be caused by tidal interactions between a massive secondary and a coplanar primary (Lubow 
&; Pringle 1993). A warping instability can be caused by irradiation from the primary (Pringlc 1996, 1997). 
A warping can be caused by direct tidal forces from a secondary that is on an inclined orbit. A disc warp 
can be caused by misalignments of the spin axis of a compact and/or magnetized primary and the disc axis. 
Of these possibilities, Murray & Armitage (1998) find that CV DN accretion discs do not seem to tilt out 
of the orbital plane by instabilities. 

If a warp is present, the warp may or may not propagate. As reviewed in Pringle (1999), an inviscid, 
locally isothermal or locally polytropic fluid disc that has Keplerian rotation, that is within a spherically sym- 
metric gravitational potential, and that is wholly tilted by the same finite angle has azimuthal wavenumber 
m=l[and thus the fluid disc has (m=l)-fold symmetry in the azimuthal direction] and radial mode number 
n=0 (and thus the fluid disc has n=0 local minima). For this disc, the wave mode frequency is zero. If, 
however, not all of the inviscid disc is tilted by the same amount, then a wave can radially propagate (see 
Papaloizou Sz Lin, 1995 and e.g., Lubow, Ogilvie, & Pringle, 2002). If the disc has some viscosity, then a 
disc wave can only propagate if the viscosity is low and the tilt angle is approximately less than H/r where 
H is vertical disc thickness and r is disc radius. If the disc has low viscosity, then the wave is transferred 
through the disc by wave action. Likewise, if the disc has high viscosity, then the wave is transferred through 
the disc by dispersion (Papaloizou & Pringle, 1983). If a quiescent, viscous disc is wholly tilted by the same 
amount and the tilt angle is greater than H/r, a wave should not propagate. 

Regardless of how an accretion disc tilts or warps and whether the wave propagates, non-coplanar 
accretion discs in binary systems precess in the retrograde direction (i.e., in the direction opposite to spin 
and orbital motion). Retrograde precessional activity has been suggested for X-ray binaries (e.g., Katz 
1973, Larwood 1998, Wijers &; Pringle 1999), for common systems that produce jets (e.g., Livio 1999), 
for protostars (e.g., Papaloizou & Terquem 1995, Larwood et. al. 1996, Larwood 1997), for CV systems 
(Bisikalo et al. 2003, 2004) and for quasars and supermassive black holes (e.g., Romero et al. 2000, Caproni 
&; Abraham 2002), to name a few. In many of these binary systems, the source to the retrograde precession 
is often suggested to be due to direct tidal forcing. In Montgomery (2009a), we show that tidal torques 
by the secondary on a tilted, spinning accretion disc is a source to retrograde precession. A more familiar 
example is the tidal torque by the Sun and Moon on the oblate, spinning, tilted Earth that causes the 
Earth's First Point of Aries to precess retrogradely. 

Of note is that some find retrograde precession occurring within a disc even though the entire accretion 
disc remains coplanar with the orbit: In numerical simulations by Bisikalo et al. (2003, 2004), a coplanar 
inner spiral density wave is generated within the inner annuli of a coplanar accretion disc and only the 
coplanar spiral density wave retrogradely precesses. Based on their numerical simulations, the wave is 



generated by sheared elliptical orbits and these sheared elliptical orbits (and thus the wave) retrogradely 
precess. Bisikalo et al. (2004) suggest that the retrograde precession is due to an influence by the mass- 
losing star by the processes discussed in Kumar (1986) and Warner (1995, 2003). However, Kumar (1986) 
finds the retrograde precession on a tilted disc's outer annulus due to the tidal influence of the secondary. 
Since the spiral density wave is found in inner annuli of the accretion disc, theory is needed on retrograde 
precession for a wave in this location. 

In Montgomery (2009a), we suggest that care must be taken in choosing the correct state of matter 
and the correct geometrical shape of the accretion disc as different shapes can produce different retrograde 
precessional values. For example, Montgomery (2009a) shows that retrograde precession of a non-magnetic, 
tilted CV DN gaseous accretion disc can best be described by a retrogradely precessing and differentially 
rotating tilted accretion disc, and superimposed on the disc is a rotating and precessing ring that is located 
near the innermost annuli of the disc. The superimposed ring is likely caused by gas stream overflow and 
additional particle migration into the innermost annuli of the tilted disc (Montgomery 2009a, b). Instead of 
a precessing inner ring in accretion discs, some accretion disc simulations show an inner spiral density wave 
in the same location: For semidetached binaries, Bisikalo et al. (2003, 2004) show a one-armed spiral shock 
wave is generated in inner annuli of cooler discs (i.e., T W^K) of binaries with high secondary-to-primary 
mass ratios. Therefore, a study is needed on whether such a wave can be generated in simulations of lower 
mass ratio, and can such a wave explain retrograde precession in CV DN systems. 

In this work, wc generate theory on retrograde precession of an accretion disc that contains an inner 
spiral density wave. We assume the wave is within the plane of a quiescent disc, yet the disc is tilted five 
degrees relative to the orbital plane. We wish to establish whether replacing the precessing inner ring found 
in Montgomery (2009a) with an inner spiral density wave has similar effects on disc retrograde precession. 
To justify the disc geometry (i.e.. disc with a coplanar inner spiral density wave), we generate a low mass 
ratio numerical simulation using the grid code of Bisikalo et al. (2003). In §2 of this paper, we develop 
the theoretical expressions to describe a retrogradely precessing, tilted CV DN accretion disc that has an 
inner spiral density wave in the plane of the disc. In §3, we describe the two numerical codes that generate 
retrograde precession in CV DN accretion discs and we present numerical simulation results. In §4, we 
compare our theoretical and numerical results with results obtained from observations and with previously 
obtained results. In §5, we discuss these comparisons, and in §6, we provide a summary, conclusions, and 
future work. 

2 Theoretical Expressions for Retrograde Precessing of Discs with an 
Inner Spiral Density Wave 

In the Montgomery (2009b) Smoothed Particle Hydrodynamic (SPH) numerical simulations of a differentially 
rotating, retrogradely precessing, tilted accretion disc, a ring forms near the innermost annuli of the disc due 
to particle migration into this inner annuli and due to gas stream overflow of the tilted disc edge. Instead 
of a ring, numerical simulations by Bisikalo et al. (2003. 2004) show gas stream build-up in the form of a 
one-armed spiral density wave. Therefore, we consider a hypothetical disc geometry for non-magnetic CV 
DN system that exhibits negative superhumps in their light curves. The hypothetical disc geometry is a 
combination of the disc geometries in both numerical simulations: We assume the disc is wholly tilted five 
or more degrees degrees relative to the orbital plane. Instead of a denser ring of material in the innermost 
annuli of the disc, we assume an inner spiral density wave that remains in the plane of the tilted disc. In 
the following subsections, we find the disc's net retrograde precessional value. 



2.1 Retrograde Precession, Dynamical Term 



The retrograde precessional angular frequency of a differentially rotating, circular, tilted disc is found to be 
(see e.g., Montgomery, 2009a) 

In this equation, q = M2/M1 is the mass ratio of the secondary-to-primary masses, 9 is the obliquity angle, 
rd is the radius of the disc, d is the orbital separation, and Uorb is the orbital frequency. This equation applies 
to a thin disc that lacks an inner hole or has a small inner hole (e.g., non-magnetic, CV DN systems). This 
disc also does not have any surface features such as spiral density waves or dense rings (or winds, jets, 
planets, holes, gaps, or any other feature that may change the angular momentum, and thus moment of 
inertia, of the system). The shape of the disc is ideally circular although the equation is assumed applicable 
to non-circular discs as well (see Montgomery 2009a). Equation (1) is valid so long as the disc processes as a 
unit and the distance between the stars d is much greater than the radius of the primary (see Montgomery 
2009a). 

The outer radius of the approximately circular disc is assumed to be 

rd = 0.6d/{l+q) (2) 

as given by Paczynski (1977). This equation is valid for mass ratios in the range 0.03 < q < 1. 

As retrograde precession is associated with tilted discs and as negative superhumps are generated when 
disc tilt is high enough (see Montgomery 2009b), then negative superhumps could be an observational beacon 
that indicates that the disc may be tilted. As the negative superhump period P— and the orbital period Porb 
are fairly precise clocks, we would like to relate the retrograde precessional period to these precise clocks. 
From observations, wc know that the orbital period is slightly longer than the negative superhump period. 
Both periods are related to the net retrograde precession period by 

If we assume a primary mass Mi=0.8Mq and mass ratio q = 0.4, then the secondary mass M2=0.32Mq, 
~ 0.43d where d ~ 1.23i?0 (Montgomery 2009b). From Newton's version of Kepler's third law, the 

orbital period and the rotation period at the disc edge is Porb = 3.852 hours and one hour, respectively. 
From Montgomery (2009b), the negative superhump period for this mass ratio is P_ ~ 3.662 hours and 
therefore, from Equation (3), Pr,net ~19.5Por6 hours, assuming small angle approximation for disc 
tilt. Notice that the spin rate at the edge of the disc is fast compared to the retrograde precessional rate. 
However, this type of precession is not fast precession because the wobble rate is not comparable to the spin 
rate. Instead, this kind of precession is driven by torques (see e.g., Montgomery 2009a). 



2.2 Retrograde Precession, Pressure Term 

To find the effects an inner spiral density wave has on retrograde precession, we first review the theory 
generated for CV DN accretion discs that progradely precess, that have spiral density waves, and that 
generate positive superhumps. From the location of the spiral density waves in a progradely precessing disc, 
we can find the location of the inner spiral density wave. 



2.2.1 Review: Net Prograde Precession in CV DN 

Spiral density waves contribute to retrograde precession through pressure as Lubow (1992) finds for pro- 
gradely precessing CV DN accretion discs that show positive superhumps in their light curves. Lubow (1992) 
finds that the net prograde precession oop-net is the net of a dynamical component, a pressure component, and 
a transient component. Of the three, the dynamical component contributes most to prograde precession. 



The second of the three is a retrograde component and is a correction term to the dynamical component. 
This pressure component is necessary to more correctly explain the net prograde precession in steady state 
discs. The third component is also a correction term and can contribute to the prograde dynamical compo- 
nent or to the retrograde pressure component. However, this transient component is most effective during 
the development stage of positive superhumps. If only quasi-steady state discs are considered, then this 
third component can be ignored. 

Pressure effects within progradely precessing discs could be due to a two-armed spiral mode. The source 
to this spiral shock generation is the j : j — 1 = 3 : 2 resonance (Lubow 1991a, b), using the terminology 
established by Frank, King, & Raine (1992). Dissipation from these spiral shocks are created from converging 
gas flows, and this dissipation is the source to the observed positive superhump as suggested by Smith et 
al. (2007). The location of these spiral density waves is in the middle and outer annuli of the quasi-steady 
state disc (see e.g.. Smith et al., 2007). 

Assuming this disc geometry for prograde precession, an analytical expression for the pressure component 
is developed in Lubow (1992). The standard Wentzel-Kramers-Brillouin (WKB) relation for spiral density 
waves in the tight winding limit is given by 



^press,outer ) 



(4) 



where Vl' and k' arc the particle angular and radial frequencies, respectively, ujpress,outer is the oscillation 
frequency of the spiral density waves, c is the gas sound speed within the disc, and k' is the radial wavenumber 
of the spiral arms (see e.g., Goldreich k, Tremaine 1979). In a non-precessing disc, fi' = n' , and in a slowly 
precessing disc, k' . In the limit where C0press,outer ^ ^' or k' and applying the previous approximation, 
Lubow (1992) approximates press. outer in the previous equation. Then Lubow (1992) couples ojpress, outer 
with the dynamical motion within the disc to find the net prograde precession: 
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(6) 



We identify the last term in the equation as the pressure term associated with these spiral density waves 
located in the middle-to-outer portion of the disc. Note that we assume quasi-static discs; the net prograde 
precession has only a dynamical prograde component and a retrograde pressure component. 
In e.g., Montgomery (2001), we show that the positive superhump period excess is 
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where we have substituted 



Pp,net — Pari i^^) 

into Equation (7) to obtain Equation (8). In these equations, P+ is the positive superhump period and Pp^net 
is the net prograde precessional period. If we assume ujp^net ^ ^orb: then we obtain the positive superhump 
period excess 
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(12) 



The pitch angle i of the spiral wave is related to the wavenumber by 



coti' = k'r' (13) 

where r' is the radial distance to that fluid element within a spiral arm. 

In Montgomery (2001), we assume average values for c and i' based on our numerical simulations. We 
find the average distance to the two numerically generated spiral density arms and the average pitch angle 
to be r' ~ 0A77d and i' ~ 17°, respectively, where d is the binary separation. After substitution of r', i', k', 
il' = 3wor&) and the Smith & Dhillon (1998) secondary mass-period relation into Equations (5) and (6), we 
find 
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2.2.2 Retrograde Precession Effects from Inner Spiral Density Wave 

As reviewed in Kato (2001), local analytical studies of hydrodynamic accretion discs have found three types 
of trapped global modes. One of the three types is inertial acoustic p-modes (the other two are gravity g- 
modes and vertical c-modes and are not discussed in this work). For a two-dimensional disc with no vertical 
modes, the p-mode is described by 

(tTiO ^press , inner) — 1^ '\' k C (15) 

where and k" are the particle angular and radial frequencies, respectively, ujpress,inner is the oscillation 
frequency of an inner spiral density wave, c is the gas sound speed within the disc, k" is the radial wavenumber 
of the spiral arms, and m is an integer that indicates the different azimuthal modes. Notice that this equation 
is very similar to Equation (4), however this equation is evaluated at the radius r" of inner spiral density 
wave and hence the wave number and particle angular and radial frequencies are different. Also note that 
we assume the same sound speed throughout the disc. 

Global spirals in discs are the m = 1 p-modes. For a one-armed spiral wave, we can reduce the previous 
equation to 

^pressAnner — ^ i \ ri^-\-k^^ . (f ^) 

In a non-precessing disc, fi" = k" , and in a slowly precessing disc, 0" k". If the disk is also geometrically 
thin and ck" <^ Q," (i.e., we have a global oscillation), then the previous equation reduces to 

j^'/2g2 

^press,inner ~ 2^}/^ ("'"'^) 

This result is similar to that found by Kato (1983). As Kato (1983) noted, this oscillation at r" moves much 
slower than the angular velocity of the disc rotation. Therefore, the oscillation acts to slow the motion of 
the gas particles in the disc and hence the sign is negative in this equation. 

To find the net retrograde precession in the disc, we need to sum the effects due to the inner spiral 
density wave with the effects due to the dynamical motions in the smooth, featureless disc: 

Wr-,net - i^r,dyn " 2^1/ V-*""/ 
~ ^r,dyn ^press,inner {^^) 

The last term in the equation is identified as the pressure term associated with the inner spiral density wave 
number k" and particle angular frequency ft". Notice that these equations are very similar to Equation (5) 
and (6), however these equations are evaluated at at the radius of the inner spiral density wave r". 



In Montgomery (2009b), we show that the negative super hump period excess is 
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Pr,net ~^ Porb 

where we have substituted Equation (3) into Equation (20) to obtain Equation (21). In these equations, 
P_ is the negative supcrhump period and Pr,net is the net retrograde precessional period. If we assume 
^p,net ^orb-: then we obtain the net negative super hump period excess 

^r.net 



e-,net ~ — (23) 
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As 0Jr,dyn is negative, then Ur^net and e-^net are also negative. If we assume a secondary mass-period 
relationship by Smith & Dhillon (1998), then we can find the net negative super hump period excess in terms 
of orbital period (not shown). 

The pitch angle of the inner spiral density wave i" is related to wavenumber k" by 
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Upon substitution into Equations (18) and (19), 
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where in the second version of ojpress,inner we have assumed r"~0.16d for the average outer radius of the 
inner spiral density wave. This location is the average outer radius of the inner ring seen in the q = 0.4 
simulation of Montgomery (2009a) and is larger than 

fcirc ~0.13d which is the minimum radius of the ring's 

outer radius. 

If d is obtained from Kepler's Third Law, c ~20 kms~^ and i" ~ 17° (Montgomery 2001), then only 
the particle orbital frequency fi" is unknown. However, we can easily find the particle orbital frequency 
through Kepler's Laws: If we solve Newton's version of Kepler's third law for the orbital period and ratio 
it with Newton's version of Kepler's third law for the period of a particle within the disc P at any radius r, 
ignoring perturbation effects by the secondary, then we find 
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After substitution for r" and q, we find P" ~ 0.08Por6 or Q," = 12. Sworfe- We can use the Smith & Dhillon 
(1998) secondary mass-period relation to find 

1-4^-1/1.58 



UJorb = 2.54 X 10 

where ^orb is in units of rad s~^. Upon substitution into Equations (26) or (27), we find 
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and C0press,inner ~ 2.1 {ojpress, outer)- As expected from Keplerian motion, the inner spiral density wave angular 
frequency is significantly different from the outer spiral density wave angular frequency. 



3 Numerical Simulation Codes and Results 

In this section, we describe two different numerical codes that generate accretion discs like those expected 
for non-magnetic CV DN. In §3.1, we summarize a grid code that generates accretion discs with an inner 
spiral density wave as shown in e.g., Bisikalo et al. (2003, 2004). In this work, we present the q = 0.4 
numerical simulation using this code to see if lower mass ratio systems can also generate an inner spiral 
density wave. In §3.2, we summarize the SPH code that generates accretion discs with a denser inner ring 
if the disc is tilted at least a few degrees as shown in Montgomery (2009a). 

3.1 3D, High Resolution, Gas Dynamical Grid Simulations 

The details of the three-dimensional, high resolution gas dynamical grid code can be found in Bisikalo et 
al. (2003). Here we provide a summary. For this work, we consider a semi-detached binary consisting of a 
donor secondary star with mass M2 = O.32M0 that fills its Roche lobe. An accretor primary star with mass 
Ml = O.8M0 is separated from the donor secondary by a distance d = 1.23i?0. The modeling is carried 
out in a non-inertial reference frame rotating with the binary and on a regular three dimensional grid in 
Cartesian coordinates. The origin of the coordinate frame coincides with the center of mass of the donor 
secondary. 

Since the problem is symmetrical about the equatorial plane, the modeling is carried out in half the 
space: A symmetrical boundary condition is imposed at the corresponding boundary of the computational 
domain. Steady-state boundary conditions are specified at the remaining boundaries of the domain: The 
density is pb = 10~^pLi, where pLi =1.5 x 10"'' g • cm~^ is the matter density at the inner Lagrange point 
Li; the equilibrium temperature of the optically thin gas is T = 13, 600^ (see Bisikalo et al. 2003); and the 
velocity is Vb = 0. The white dwarf accretor primary is specified to be a sphere with radius Ri = O.OlSd, 
and a free-inflow condition is imposed at its surface: All the matter that falls in cells occupied by the white 
dwarf accretor is considered to have fallen onto the star. 

The force field in the system is described by the Roche potential: 

GM2 GMi 
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(32) 



where G is the gravitational constant and Q is the rotational speed of the system. The gradient of the Roche 
potential is equal to zero at five libration points (the so-called Lagrange points): V<I> = 0. The equipotential 
surface that passes through the inner Lagrange point Li forms the Roche lobes of the system's components. 

The size of the donor secondary, a red dwarf, is specified by the boundaries of its Roche lobe. The 
following boundary conditions have been selected at its surface: p2 = pli, ^2 = 3200 K, V2 = vli • V$/|V<I>|. 
The speed at the inner Lagrange point Li is specified to be the local sonic speed: = Cs=6.6 kms~^. 

The region within the donor secondary is excluded from the calculations. 

We describe the flow structure in the binary using a system of gravitational gas dynamic equations, 
taking into account the radiative heating and cooling of the gas in the optically thin case: 

f |e + v.(H = 0, 

^ + V-{pv0v) + VP = -pV^-2[nxv]p, ^^^^ 
M^WM. + V • {pv){e + P/p + z;V2) = 

- pvv^ + p^mf [r(r, T„) - A(r)] . 



Figure 1: The inner spiral density wave is shown in cartesian coordinates as fractions of the separation 
distance d for the q=0.4 numerical simulation. Note the Z axis scaling is significantly smaller than the 
others, and only the upper half of the disc is shown. The secondary is not shown. However, the secondary 
is orbiting counterclockwise, and it is located at (X=0, Y=0, Z=0). 

Figure 2: The inner spiral density wave from Figure 1 is shown for orbits 30, 30.27, 30.49, and 30.74. 
To convert into time, multiply one orbit by Port- Cartesian coordinates are in fractions of the separation 
distance d. The secondary is orbiting counterclockwise and is located at (X=0, Y=0, Z=0) in all panels. 

Here, p is density, v= {u,v,w) is velocity, P is pressure, t is time, e the inner energy, rup is the mass of the 
proton, and T{T,Ta) and A(T) are the radiative heating and cooling functions (see Bisikalo et al. 2003). 
The system of the gasdynamic equations is closed with the equation of state for an ideal gas, P = (7 — l)pe, 
with an adiabatic index 7 = 5/3 for a monatomic ideal gas. 

We solve the system of equations using the Roe-Osher-Einfieldt method adapted for multiprocessor 
computers (Roe 1986; Chakravarthy & Osher 1985; Einfeldt 1988). This technique makes it possible to 
obtain high-accuracy solutions of gas dynamical problems. The simulations are carried out until the system 
is in the steady-state accretion regime. 

Figure 1 shows that a spiral arm can be generated in the steady state solution for the lower mass ratio 
q = 0.4 numerical simulation where Mi = 0.8Mq. The cartesian coordinates (X,Y,Z) are in fractions of 
the separation distance d. Although in this work we only show half the space, the full three-dimensional 
simulations also have the same spiral wave as shown in the purely hydrodynamic solutions by Bisikalo et 
al. (2008) and in the hydrodynamic solutions with magnetic fields by Zhilkin et al. (2010). Therefore, the 
spiral density wave is not a consequence of symmetry and not due to magnetic fields. 

As shown in Figure 1, the spiral density wave does have a small projection relative the disc mid-plane. 
Figure 2 shows snapshots of the disc over three-quarters of an orbit. In their numerical simulations, Bisikalo 
et al. (e.g., 2003, 2004) find the wave to be a result of sheared elliptical particle orbits in the innermost 
annuli of the disc: As the particles travel in elliptical orbits around the primary, where the primary is at 
one common focus to the elliptical orbits, the apastron points of the elliptical orbits are gravitationally 
torqued by the secondary to precess in the retrograde direction. The closer an apastron point is to the 
secondary, the larger its semi-major axis, and the more torque that apastron point experiences. Because 
velocity flow is very slow at these apastron points and because flow lines of the gas stream cannot cross, 
gas particles accumulate near these apastron points. We find in this work that these accumulations have 
two sources: gas particle migration from outer annuli and, as shown in Figure 10 of Bisikalo et al. (2005), 
the gas stream that overflows the disc edge and travels radially across the disc face to inner annuli. The 
result of an over population of gas particles near apastron points is an elliptical spiral density wave in the 
innermost annuli. Because apastron points are not uniformly spaced in angle due to differential rotation of 
the disc, accumulation of gas particles is not uniform near these points. The result is a variation in density 
along the wave. The net of the two effects is an asymmetric, elliptical, tilted spiral density wave located in 
inner annuli as shown in Figure 1. Because the center of mass is not in the center of the elliptical spiral and 
not in the orbital plane, the spiral density wave is subject to torques by the secondary and thus the spiral 
density wave retrogradely precesses. We find for the q=0.4 numerical simulations in steady state that as 
the secondary orbits the center of mass, the secondary returns to its starting position (relative to the spiral 
density wave) at a period that is ~3% shorter than the orbital period. This shorter period indicates that 
the spiral density wave is in retrograde precession. In ~ 33 orbits, the spiral density wave has made one full 
retrograde turn. We note that the disc does not retrogradely precess. 

From Figures 1 and 2, we find that the inner spiral density wave outer radius ranges from ~ (0.1 — 0.2)d. 
However, the average outer radius is ~0.16d, confirming our assumption that the inner spiral density wave 
is located near the inner ring seen in the SPH simulations Montgomery (2009a). 



3.2 3D, 100000 particle, SPH Simulations 



For the numerical simulations in Montgomery (2009b), we use Smoothed Particle Hydrodynamics (SPH) 
which is a Lagrangian method that models highly dynamical, astrophysical fluid flow as a set of interacting 
particles (see Monaghan 1992 for a review). The code utilizes the simplest SPH form - constant and uniform 
smoothing length h and constant and uniform mass particles. Modifications to the original code include 
increasing particle number to 100,000 and adopting a Smith & Dhillon (1998) secondary-to-primary mass- 
period relation. The code is discussed further in Montgomery (2009b), however applicable information is 
also summarized here. 

The original code models hydrodynamics assuming an ideal gamma-law equation of state P = (7 — l)pu 
where P is pressure, 7 is the adiabatic index, p is density, and u is specific internal energy. The sound speed 
is Cs = y/ji'J — 1)'W. The momentum and internal energy equations per unit mass are, respectively, 

£v VP^„ GMi GM2 

^ = -— + f-c - -^r, - -3-r2, (34) 

and 

du P„ 

-77 = V • V + e^isc (35) 

at p 

in their most general form (Simpson 1995). In these equations, fyj^c is the viscous force, e^^sc is the energy 
generation due to viscous dissipation, and ri and r2 are the displacements from stellar masses Mi and M2, 
respectively, that are given in units of solar masses. The momentum and energy equations in SPH form for 
particles i and j are, respectively. 



'i2 

and 

dui 



GMi GM2 

— 3— ra 3— rj2 (36) 



dt ■ ('^^ 



or, if the last equation fails in preventing negative internal energies, 

dui Pi 



In these equations, Hjj is the Lattanzio et al. (1986) artificial viscosity 



^Y.mj{l + Uij)vij-ViWij. (38) 



^ I -apij + Pfil Vij ■ r,.j < 

'•^ I otherwise 



where 



= ^^2^?2^ (40) 

and Cg^ij = ^{cs,i + Cgj) is the average sound speed, Vjj = Vi — vj, Vij = Vi — rj, and r/^ = O.Ol/i^ as shown 
in Simpson & Wood (1998). 

For our numerical simulations, we choose a = ^ = 0.5. Our viscosity is approximately equivalent to a 
Shakura & Sunyaev (1973) viscosity parametrisation (1/ = a'cgH where H is the disc scale height) a'=0.05. 



As Smak (1999) estimates a' ~ 0.1 - 0.2 for DN systems in high viscosity states, our simulations are more for 
quiescent systems. As only approaching particles feel the viscous force and since neither radiative transfer 
nor magnetic fields are included in this code, all energy dissipated by the artificial viscosity is transferred 
into changing the internal energies. As radiative cooling is not included in this code, an adiabatic index 
7=1.01 is incorporated to prevent the internal energies from becoming too large. 

By integrating the changes in the internal energies of all the particles over a specific time interval n, 
variations in the bolometric luminosity yield an approximate and artificially generated light curve 

L„ = ^d<. (41) 

i 

All simulations start by injecting five particles per major time step, 200 major time steps per orbit, or a total 
of 1000 particles per orbit. The injection velocity is based on cin injection temperature T^fij — 4x10"^ K and 
is approximately the sound speed that has been scaled to system units. The code adopts an approximate 
secondary mass-radius relation R2 = {M2/Mq)-'^^I'^^Rq that applies for 0.08 < < 1.0 (Warner 2003). 
In this relation, R2 and M2 are in solar radii and solar mass, respectively. 

A unique feature of this code is the conservation of disc particle number. Any time a particle is accreted 
onto the secondary or primary mass or is lost from the system, a new particle is injected through the 
inner Lagrange point Li. We build a disc of 100,000 particles and maintain this number throughout the 
simulation. After building the disc, we then allow the disc to evolve in the short term to a quasi-equilibrium 
state where only a few particles are injected per orbit. At orbit 200, we artificially rotate the disc out of the 
orbital plane 5°. 

For our simulations, we assume a primary mass Mi=O.8M0 and we vary the secondary mass M2 such 
that 0.35 < q < 0.55. For example, ii q = 0.4 then M2=0.32 Mq and our simulation unit length scales 
to d ^ 1.23i?0 where we have used the secondary mass-radius relationship shown above and the Eggleton 
(1983) volume radius of the Roche lobe secondary 

^ = (42) 

d 0.6g2/3 + ln(l + gV3) • ^ ' 

The Eggleton (1983) relation is good for all mass ratios, accurate to better than 1%. Taking the radius of 
the primary to be 6.9x10^ cm, the scaled radius of the white dwarf is i?i = 0.0081ci. 

After building the disc, we allow the disc to evolve in the short term to a quasi-equilibrium state where 
particles are injected at the rate they are removed from the system by either being accreted onto the primary 
or the secondary or lost from the system. The average net rate of accretion for a quasi-static disc is around 
500 particles per orbit as shown in Montgomery (2009b). As this net rate is approximately half that injected 
to build the disc, the steady state mass transfer rate reduces to m 1 x 10~^^ Mq/jt or a rate similar to 
an SU UMa in quiescence. We estimate the density of the gas near Li to be p ~ 10~^° g cm~^ using 
nip = (4/3)7r/i^/9. As we maintain 100,000 particles in the steady state disc, then the total mass of the disc 
is Maisc= 100,000 x nip = 3.5x10^^ g or M^isc ^ 2xlO~^^M0, a negligible value compared to the mass of 
either star. If we assume the outer radius of the disc is Rdisc ~ ^fdrc, or twice the circularization radius 
(Warner 2003) 

-0.426 



^ = 0.0859g-"-^"'^ (43) 
a 

where 0.05 < q < I, then for q = 0.4 we find rdrc ~ 0.13d and Rdisc ~ 0.25d. As p « MdiscRd^sc ^^^^ 
gravity is neglected in these simulations. 

As output to the code, we create an artificial light curve for each simulation (see Montgomery, 2009b). 
A non-linear least squares fit to each Fourier transform is computed and the standard deviation to twice 
the negative superhump frequency 2z^_ is found. This data in Table 1 is from Montgomery (2009b) and, 
along with the calculated negative superhump period, the propagated error and the negative superhump 
period excess are also listed. The units for period are hours whereas the units for frequency are orbits"^ to 



Table 1: Negative Super hump Simulation Data 



Frequency (Pgrb) Period (hr) 



q 


2i/_ 




Porb (lir) 


P_ 


0"P_ 




0.35 


2.098 


0.004 


3.540 


3.374 


0.006 


0.0467 


0.375 


2.105 


0.004 


3.698 


3.514 


0.006 


0.0498 


0.4 


2.104 


0.004 


3.852 


3.662 


0.006 


0.0494 


0.45 


2.110 


0.004 


4.150 


3.934 


0.007 


0.0521 


0.5 


2.115 


0.003 


4.436 


4.195 


0.007 


0.0544 


0.55 


2.119 


0.003 


4.712 


4.447 


0.007 


0.0562 



be consistent with the Fourier transform figure. To convert to real units, multiply one orbit by the orbital 
period. 

3.3 Grid Versus SPH Numerical Methods 

Both the SPH simulations and the grid simulations have gas particle migration from outer annuli to inner 
annuli of the disc and have gas stream overflow. However, the SPH simulations have an imposed disc tilt 
whereas the grid simulations do not. The SPH simulations generate an inner ring of higher density whereas 
the grid simulations generate an inner spiral density wave. The SPH simulations have gas particles traveling 
in circular orbits in inner annuli whereas the grid simulations have gas particles traveling in elliptical orbits 
in inner annuli. Thus a wave is produced in the grid simulations but not the SPH simulations. 

SPH is a viscous method that heavily relies upon particle interaction. If particle neighbors are few then 
fine resolution is lost. Grid simulations are also not immune to resolution: To see detail like an inner spiral 
density wave, a fine grid is needed like the one used in this work. SPH also heavily relies on artificial viscosity. 
The smoothing effect of the artificial viscosity is known to hinder the production of spiral waves (see e.g., 
Lanzafame 2003). Although both the total particle number and the artificial viscosity may contribute to 
particle orbit shape in inner annuli of SPH simulations, the latter is likely the leading cause to the lack of 
an inner spiral wave in the SPH simulations. 

4 COMPARISONS WITH OBSERVATIONS 

Table 2 lists observational data including orbital periods; negative and positive superhump periods; and 
negative and positive superhump period excesses for several CV systems. If errors are known, they are 
listed in parentheses. We select long orbital period systems so that we may compare their observational 
results with our higher mass ratio numerical simulation and with our theoretical results. In the table, NL 
is novalike, PS is permanent superhumper, and IP is intermediate polar. 

The disc geometry in this study is hypothetical as it is a combination of the disc geometries found in 
numerical simulations of Montgomery (2009b), Bisikalo et al. (2003, 2004), and this work. The numerical 
simulations of Montgomery (2009b) show a tilted disc with a coplanar denser inner ring. The numerical 
simulations of Bisikalo et al. (2003, 2004) and this work for q = 0.4 show a coplanar disc with a dense 
inner spiral density wave that projects at a slight angle to the disc mid-plane. Hence the hypothetical disc 
geometry for this study is a fully tilted disc around the line of nodes and the disc contains an inner spiral 
density wave that remains coplanar with the disc. 

Figure 3 shows a comparison of the numerical negative superhiimp period excess (A) to the observational 
negative superhump period excess (*) values listed in Tables 1 and 2, respectively. In this figure, the width 
of the symbol can be taken as the error in the orbital period and, to express mass ratio in terms of orbital 
period, we include the Smith &; Dhillon (1998) secondary mass-period relation. As the geometry of the 
accretion discs for each observational system is unknown, we do not differentiate between different CV types 
with different symbols even though IPs probably have a significant absence of inner disc annuli, NLs could 



Tabic 2: Positive and Negative Superhump Observational Data, Calculated Precessional Periods, and Period Excesses 
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Figure 3: Negative Superhump period excess as a function of orbital period. Observational data are shown 
as stars and numerical data are shown as triangles. The dashed line represents dynamical only contributions 
to retrograde precession. The dotted line represents a spinning, tilted disc with an inner spiral density wave. 
The short dash line combines the effects of a differentially rotating, tilted disc with the pressure effects from 
an spiral density wave. The solid line combines the effects of a differentially rotating, tilted disc and an 
attached rotating innermost ring. Data are from Tables 1 and 2. disc is tilted 5 degrees. 

have discs with spiral density waves that cause the positive superhump as well as an inner spiral density 
wave or ring that causes the negative superhump, SU UMa's could have only the spiral density waves that 
cause the positive superhumps, etc. Therefore, caution is advised when comparing the observational data 
with the numerical simulations as the scatter in the observational data is likely due to different geometries 
of a fluid disc (see Montgomery 2009a). We also note that we have not taken into account other factors 
that may affect angular momentum and retrograde precession such as mass loss from winds, jets, mass loss 
to circumbinary rings, etc. 

In Figure 3, we plot as a dash-dot line a differentially rotating, retrogradcly precessing disc [i.e., from 
Equation (1)]. The dashed line represents a differentially rotating ring that can be found near the innermost 
annuli of the disc as discussed in Montgomery (2009a). The dotted line represents Equation (31), the 
retrograde precession due to an inner spiral density wave that is located near the ring's radius. As shown, 
the dotted line compares well with the dashed line. Likewise, the net effect of a differentially rotating disc 
with a spiral density wave (the short dash line) compares well with the net effect of a differentially rotating 
disc with an attached, inner ring (the solid line). In other words, the effects of pressure within the disc due 
to an innermost spiral density wave or the effects due to an attached innermost ring are necessary to better 
explain net retrograde precession in CV DN tilted, accretion discs should either be present in inner annuli. 

For the q=0.4 numerical simulation presented in this work, we find that the spiral density wave takes 
~33 orbits or '^5.297d to retrogradely precess. The negative superhump period is ~3% shorter than the 
orbital period or ~3.736 hours. The negative superhump period excess is ~0.03. 

5 DISCUSSION 

Prom the results, we can say that non-magnetic CV DN that show negative superhumps in their light curves 
can be described by a tilted, differentially rotating, precessing disc (i.e., the dynamical component) that 
has either a retrogradcly precessing spiral density wave located in the innermost annuli or a retrogradely 
precessing ring that is superimposed on the innermost annuli of the disc. Therefore, with the addition of 
either a retrogradely precessing ring or a retrogradely precessing spiral density wave, the geometry of the 
disc has changed as well as the net angular momentum and the net retrograde precession due to the change 
in the moment of inertia. Unknown is whether an inner spiral density wave or an inner denser ring is present. 
However, both the ring and the spiral density wave have the same effect on the disc so long as both have 
nearly the same radial extent and are coplanar with the tilted disc. If either is not coplanar with the disc, 
then differential precession would be present, effects not seen in non-magnetic CV DN systems. 

We note in Figure 3 that we assume a Smith &: Dhillon (1998) secondary mass-period relation to obtain 
the negative superhump period excess as a function of orbital period. By assuming this relation, the secon- 
daries are taken to be main-sequence stars. Those observations that agree with our numerical simulations 
and the solid and dashed lines may have tilted accretion discs with either an innermost precessing ring as 
discussed in Montgomery (2009a) or an innermost spiral density wave (this work). 

Bisikalo et al. (e.g., 2003, 2004) find that their numerically simulated inner spiral density wave retro- 
gradely precesses yet the rest of the disc does not. We confirm this result in the q=0.4 numerical simulation 
presented in this work. Because this spiral density wave has its center of mass off-center and above the 
orbital plane of the disc, the spiral density wave is subject to torques by the secondary. The net result is 
a retrogradely precessing spiral density wave. The values for the negative superhump period, retrograde 



precession, and negative superhump period excess are long but not unreasonable when compared with the 
values in Table 2. 

6 SUMMARY AND CONCLUSIONS 

Negative superhumps are found in observational systems with high mass ratios, long orbital periods, and 
high mass transfer rates. In Montgomery (2009b), we find that the negative superhump can be due to extra 
light being emitted from inner disc annuli. The extra light is modulated by gas stream overflow that strikes 
inner disc annuli and by gas particle migration into inner annuli as the secondary orbits. For this result to 
be true, we find that the disc must be fully tilted around the line of nodes. 

In non-magnetic binary systems, discs that are tilted also retrogradely precess. Patterson et al. (1993) 
suggest that retrograde precession in CV DN accretion discs may be by the same source that causes the 
Earth to retrogradely precess. In Montgomery (2009a), we generate theoretical expressions for a disc that 
is misaligned with the orbital plane to test this suggestion, and we find that spinning, tilted, precessing 
accretion discs in non-magnetic binary systems can be described by a differentially rotating and precessing 
disc that has a precessing, denser ring superimposed on disc near the innermost annuli of the disc. In this 
work, we consider a differentially rotating and precessing accretion disc that has an inner spiral density wave 
instead of a ring but having the same radial extent as the ring and in the same plane as the disc. That is, 
in addition to dynamical effects, we consider pressure effects in our analytical expressions. We adopt the 
theory used to describe pressure effects within CV DN discs that progradely precess. 

We compare our analytical expressions generated in this work with the analytical expressions generated 
in Montgomery (2009a) and the numerical simulations generated in Montgomery (2009b) for the mass ratio 
range 0.35 < q < 0.55. In Montgomery (2009a), we find that non- magnetic, differentially rotating, tilted 
accretion discs can be described by a differentially rotating, retrogradely precessing accretion disc that has a 
superimposed retrogradely precessing ring in the innermost annuli of the disc. We find in this work that these 
same non-magnetic, differentially rotating, tilted accretion discs can equally be described by a differentially 
rotating, retrogradely precessing accretion disc that has an innermost spiral density wave so long as the 
radial extent of the wave is like the radial extent of the ring. Therefore, we conclude that pressure effects in 
retrogradely precessing systems would need to be included if negative superhumps are present in the light 
curve and if the fluid disc has an inner spiral density wave. 

As part of this work, we generate a q = 0.4 accretion disc with an inner spiral density wave using a 3D, 
high resolution, gas dynamical grid code. The inner spiral density wave is not generated due to magnetic 
fields. It is a consequence of elliptical particle orbits. The inner spiral density wave projects at a small 
angle relative to the orbital plane. The tilt is due to differing densities along the spiral density wave. As a 
result of the center of mass being off-center and not in the orbital plane, the spiral density wave retrogradely 
precesses due to tidal torques, yet the disc does not. The spiral density wave takes approximately 33 orbits 
to retrogradely precess, a value that is long but not unreasonable compared to the observations. 
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